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Abstract 

We consider the steady states of a harmonic oscillator coupled so strongly to a two-level system 
(a qubit) that the rotating wave approximation cannot be made. The Hamiltonian version of this 
model is known as the E ® (3 Jahn- Teller model. The semiclassical version of this system exhibits 
a fixed point bifurcation, which in the quantum model leads to a ground state with substantial en- 
tanglement between the oscillator and the qubit. We show that the dynamical bifurcation survives 
in a dissipative quantum description of the system, amidst an even richer bifurcation structure. 
We propose two experimental implementations of this model based on superconducting cavities: a 
parametrically driven nonlinear nanomechanical resonator coupled capacitively to a coplanar mi- 
crowave cavity and a superconducting junction in the central conductor of a coplanar waveguide. 



PACS numbers: 85.85.+j; 42.50.Wk; 82.40.Bj; 84.40.Dc; 85.25.-j 



I. INTRODUCTION. 



Circuit quantum electrodynamics [T] and quantum nanomechanicsj3] have emerged in the 
last few years as new experimental contexts for the study of strongly coupled quantum 
systems. A superconducting coplanar microwave cavity can support very strong electric 
fields with very low dissipation. The electric field due to a single photon can be as large 
as 0.2 V/m and cavity quality factors as high as 10 6 have been obtained[2j. Large electric 
dipoles, in the form of superconducting junction elements, can be placed into the gap between 
the central conductor and coupled into the cavity field which is treated as a simple harmonic 
oscillator. The coupling strength can now be made far larger than equivalent experiments 
in atomic physics, and may yet be made still larger [I]. Alternatively, a nanomechanical 
resonator can form one plate of a capacitor coupling a coherent driving field to the cavity 
field |H [5]. If the nanomechanical element exhibits a significant Duffing nonlinearity then 
under parametric driving it can be approximated as a two level system interacting very 
strongly with the cavity resonator [7J. In this paper we consider a model in which the 
coupling between the field and the two level system is so large that we cannot make the 
rotating wave approximation. The model that results is known as the E ® (3 Jahn- Teller 
model. 

The E ® (3 Jahn- Teller model describes the interaction between a single simple harmonic 
oscillator and a two level system, or qubit. It has recently received some attention as there is 
a critical coupling strength at which the nature of the ground state undergoes a morpholog- 
ical change reflecting a bifurcation in a fixed point of the corresponding classical model [6] . 
As the coupling strength increases the ground state entanglement increases monotonically. 
Starting from a zero coupling strength and increasing to one infinitely large, the ground 
states change from: 

\0) c \0) q ^\a) c \+) q + \-a) c \-) q (1) 

where |0) c and |0) are the bare cavity oscillator and qubit ground states and |±a) c is 
an oscillator coherent state while |±) are orthogonal qubit states. The rate of change of 
entanglement as a function of the coupling strength is greatest for coupling strengths near 
the critical coupling strength for a fixed point bifurcation in the corresponding semiclassical 
description [6]. This has a significant implications for the ability to reach the zero photon 
state in the cavity by cooling. If one were to engineer a system with a coupling strength 
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above the critical value, cooling would reach a ground state in which the number of photons 
in the field was \a\ not zero, while the qubit would be measured to be in a totally mixed state 
for a only a little large than unity. These statements of course apply only to Hamiltonian 
systems, without damping. As real systems have finite damping, no matter how small, it is 
our objective here to determine to what extent the ground state bifurcation exhibited in the 
conservative system is manifest in the damped system and further to specify experimental 
scenarios in circuit QED in which it may be observed. 

The paper is structured as follows. In section I we present a detailed analysis of the fixed 
point structure of the dissipative Jahn- Teller E®(3 model in a semiclassical description. We 
include both dissipation of the oscillator and the two-level system. Surprisingly, despite the 
large number of parameters in the mode, the bifurcation structure is shown to depend on 
only three dimensionless independent parameters. In section II we consider the quantum 
version of the model in which dissipation is described using a Markov master equation of 
the oscillator and the two-level system. We numerically determine the steady state of the 
system. After tracing out the two-level system, we construct the Q-fuction for the oscillator. 
We then show that as the control parameters are varied through the values at which the 
semiclassical model shows bifurcations the Q-function changes from single peaked to double 
peaked with support on the semiclassical fixed points. In section III we present two physical 
systems in circuit quantum electrodynamics and quantum nanomechanics that could be used 
to implement the dissipative Jahn- Teller model. We suggest a number of key experimental 
signatures of the bifurcation. Finally in section IV we summarise our results and suggest 
new directions for further work. 

II. THE DISSIPATIVE E ® (3 JAHN-TELLER MODEL. 

We consider the case of a two level system coupled to a simple harmonic oscillator. The 
coupling is linear in the oscillator displacement and represents a state dependent constant 
force acting on the oscillator. The two-level system Hamiltonian includes a term which mixes 
the eigenstates of the conditional displacement. In order to model a realistic device we also 
need to include dissipation of both the oscillator and the two level system. The oscillator is 
damped into a zero temperature heat bath with an amplitude decay rate of k. The two-level 
system is assumed to undergo spontaneous emission at rate 7, and phase decay at rate r. 
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The irreversible dynamics due to these processes will be described by a master equation for 
weak damping (and with the rotating wave approximation for the system-bath coupling) , 
for density matrix, p, of the total resonator and qubit system. 



^ = — ^ [H, p] + k (2a pa* — o)ap — po)a) 

7 r 



(2) 



where: 



and: 



H = hcua^a + — — a x + — a z + hX (a + a*) a z + hrj (a + a*) (3) 

= ^ i?z T i^ y ) ( 4 ) 

where a x etc. are Pauli matrices. Note that we have modelled dissipation of the qubit as 
spontaneous emission in the eigenbasis of a x . This is based on the assumption that the free 
Hamiltonian of the qubit is simply proportional to a x . This makes our model consistent with 
the Hamiltonian model discussed in [6] which has no dissipation and e — 0, and ensures that 
for this limit, the fixed points of the two models coincide. The coupling between the qubit 
and the oscillator is modelled by the term proportional to A and represents a state dependent 
force acting on the oscillator. Alternatively, we can think of this term as describing a 
dependance of the energy eigenstates of the qubit on an oscillator degree of freedom, as 
occurs in the orginal Jahn- Teller model in which the electronic energy levels are dependent 
on one or more relative nuclear coordinates. Finally we have added an independent resonant 
force acting on the oscillator degree of freedom through the term proportional to rj. In the 
circuit QED realisation this would represent a driving voltage applied to the co-planar cavity. 

In the absence of dissipation, the semiclassical equations of motion that follow for the 
Hamiltonian have a fixed point pitch-fork bifurcation [6j at a critical coupling strength of 
Kr = y ^ L - A single stable elliptic fixed point, with zero cavity field amplitude, below the 
bifurcation changes stability to give two new elliptic fixed points with equal and opposite 
cavity field amplitude. When dissipation is included we expect that a similar bifurcation will 
occur but in this case the fixed points will be zero dimensional attractors. We can study this 
bifurcation by deriving the semiclassical equations of motion as follows. Using the master 
equation we construct moment equations for the expectation of each of the two-level system 
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operators (a X) a y , and a z ) and for the quadrature phase field operators defined by 

X = Ua + ct) 

2 1 (5) 
Y = -i- (a - at) 

which satisfy the commutation relations, [X,Y] = i/2. The equations of motion for the 
expectations of the five quantities of interest are found to be: 

d<x> 



d<y> 



= u> (Y) — k (X) 
= -rj-X{a g )-u{X)-K(Y) 



dt 

^ = -e (a y ) - 4A (Xa y ) - 7 (1 + (a x )) (6) 
^ = -A (a z ) + e (a x ) + 4A (Xa x ) - (a y ) 

The equations of motion for the first order moments couple to the second order moments. 
We define the semiclassical equations by factorising these second order moments to get a 
closed system of equations. Specifically, this means that we make the two assumptions 
that (Xa x ) = (X) (a x ) and (Xa y ) = (X)(a y ), or equivalently that (X,a x ) <C (X) (a x ) 
and (X,a y ) <C (X) (cr y ). (An interesting special case arises when A = 0. Then equations 
of motion for (X), (Y), and (a z ) are seen to decouple from those for (a x ) and (a y ) and 
consequently higher order moments. This decoupled system of (X), (F), and (a z ), can be 
exactly solved). After factorising the second order moments, the semiclassical variables are 
defined by: 

(X) 



L x (7) 



(<Ty) ^ Ly 

(a z ) i ^ L z 
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Yielding the semi-classical equations of motion (where a dot indicates a time derivative): 

x = coy — kx 

y = —r\ — XL Z — ujx — Ky 



L x = -eL y - A\xL y -7(1 + L x ) (g) 

7 + r 

L y -AL Z + eL x + iXxL x — L y 

7 + r 

L z = AL y - L z 

whose steady states must satisfy the Bloch sphere constraints, which depend on the presence 
or absence of qubit dissipation. If there is no qubit dissipation, the steady states lie on the 
Bloch sphere; whereas with the presence of qubit dissipation they may lie inside: 



L 2 x + L 2 y + L 2 z = l if 7 = r = o 

L 2 X + L\ + L\ < 1 if 7 > or r > 

A. Semi-classical fixed points locations 



(9) 



The semi-classical equations of motion, ([8]), have fixed points (x = y = L x = L y = L z = 0) 
that satisfy the Bloch sphere constraints, ([9]). These solutions are markedly qualitatively 
different depending on the presence or absence of qubit dissipation. In fact, there are three 
qualitatively different semi-classical steady states: no qubit dissipation; qubit dissipation 
consisting of only dephasing (no spontaneous emission); and qubit dissipation with spon- 
taneous emission. The possible presence of oscillator/cavity decay is included in each of 
the three categories. The previous coupling positivities (cj,Z\,A > and ft, 7, r > and 
e, r] G R) will be assumed for all of the following semi-classical analysis. 

We also define several convenient parameters: first, the bias / driving parameter £, which 
we see is zero in the case of an external voltage v = — second, a parameter a dependent 
on the coupling A; third, a parameter (3 dependent on the magnitude of the qubit dissipation 
parameters; fourth, a parameter S dependent on the ratio between the two different types 
of qubit dissipation (spontaneous emission 7 and dephasing r); fifth, we also define two 
combinations of these parameters: /i and v. The previous assumptions about coupling 
positivities imply similar assumptions about these parameters {a ) \i > and /?, 5, v > 1 and 
£€R). 



4X 2 ui 



a 



/i = 5a = 



A (uj 2 + H?) ' 

4(1+5) x2 " 



13 = 1 + 



7+T 
2A 



6 = 1 + 



r 

7 




(12) 



We must also consider the stability of the fixed points. The five semi-classical variables 
can be considered as a vector x, such that about a fixed point x° we have: 

iT 



= x - x° = [x - x°, y - y°, L x - L° x , L y - L° y , L z - L° z ] 



— 5x = M<5x 

dt 



(13) 
(14) 



where the Jacobian matrix M is: 



M = 



— K 
—Ul 



— K 








-4AL° 



-7 

4AL° e + A\x° 







-e - 4Ax° 

7 +r 





-A 


-A 

i+r 
2 J 



(15) 



Stability of the fixed point requires all the eigenvalues of the Jacobian to have a real part 
less than or equal to zero. In general stability depends on many more coupling parameter 
combinations than those which define the fixed points. Thus, the stability is generally 
calculated numerically for fixed points with specific values of all couplings. 



1. No qubit dissipation 

When qubit dissipation is neglected (7 = r — 0, k > 0), there are three classes of 
semi-classical steady states. Two of these classes require £ = 0; and the other requires 

Class 1: for £ = there are two fixed points that occur for all parameter values: 



x° 


T]UJ 


e 


LU 2 + K 2 


~4A 


y° 






UJ 2 + K 2 






= ±1 






= 






= 





(16) 
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The fixed point with L° = +1 is unstable for almost all[*] coupling values; while the fixed 
point with L° = — 1 is almost always stable for a < 1 and almost always unstable for a > 1. 

Class 2: also for £ = there are two fixed points that only occur for a > 1 (note that for 
a — 1 this second class of fixed points is also the first class just described): 



o t]uj Xuj I 1 e Xuj 

UJ 2 + K? UJ 2 + K 2 V O 2 4A UJ 2 + K 2 V Of 2 



q /ay ftA / 1 fte ftA / 1 

' ! uj 2 + k 2 uj 2 + ft 2 V <^ 2 4Acj cj 2 + /^ 2 V Qf 2 



1 

a 



(17) 



l; = o 



These fixed points can be stable or unstable depending on the coupling values (for example 
they are almost always (> 99%) stable for k = 0). The L z components of the first two 
classes of fixed points are plotted in figure [T} 



[*] Stability was determined numerically. This involved taking random samples in the parameter space 
spanned by lj, Z\, e, A, 77, ft, 7, and r. "Almost always" stable/unstable indicates that the vast majority 
but not all (> 99% but < 100%) of the random samples taken yielded stable/unstable fixed points. 
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Lz(a,£=0) for no qubit dissipation 




FIG. 1: L z component of the semi-classical steady states as a function of the parameter a for no 
spontaneous emission or dephasing of the two- level system (j = T = 0) and no driving (£ = 0). 
There are two solutions along the line L? z — 0, one of which is almost always stable for a < 1; 
otherwise these L\ — solutions are almost always unstable. The two new solutions which appear 
for a > 1 have L z components L\ — ±y^l — ^ and can be stable or unstable depending on the 
coupling values (for example they are almost always stable when there is no oscillator decay). 
Hence, depending on the values of the coupling parameters there is often a supercritical pitchfork 
bifurcation at a — 1. 



Class 3: for £ ^ there are up to four real fixed points dependent on a quartic equation: 

CJ 2 + K 2 LO 2 + Hv 2 Z 



y o = ^ l q 

CJ 2 + K 2 LO 2 + Hi 2 Z 

L o = _lJL_ (18) 



X 



L° Z = L° Z 

where L° z satisfies the quartic equation (from ([9])): 

(L°) 2 + a 2 ((L°) 2 -l)(L°-e) 2 = (19) 

Note that for £ ^ 0, L° z = t; is never a solution to this equation and so the pole in the 
expression for L° above is never encountered. 

The L z component of this third class of fixed points is plotted in figure [2j The bifurcations 
of these fixed points are shown in the bifurcation diagram of figure [3| 
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Lz(ff,^^0) for no qubit dissipation 




FIG. 2: L z component of the semi-classical steady states as a function of the parameters a and 
£ for no spontaneous emission or dephasing of the two-level system (7 = r = 0) with driving 
(£ 7^ 0). There are up to four solutions where the L z components are the real roots L\ of the 
quartic equation (L®) 2 + a 2 ^(L^) 2 — 1^ [LP Z — £) 2 = 0. It is clear that varying either aor( can 
take the solutions through bifurcations. This is shown explicitly in the bifurcation diagram of 
figure |3j 



2. Dephasing- only qubit dissipation 

When the qubit dissipation is considered to consist of only phase decay (7 = 0, r > 0, 
k > 0), there are two classes of semi-classical steady states. One of these classes requires 
£ = 0; and the other requires £ 7^ 0. 

Class 1: for £ = there are infinite fixed points that occur for all parameter values: 



x° 


T]UJ 


e 


UJ 2 + K 2 




y° 


KT] 


K€ 


UJ 2 + K 2 


AXuj 


L° x 


= L° X 






= 






= 





where (LP X ) < 1. These fixed points can be stable or unstable depending on the coupling 
values and the choice of L x . 
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FIG. 3: Bifurcation diagram of the semi-classical steady states as a function of the 
parameters a and £ for no spontaneous emission or dephasing of the two-level system 
(7 = r = 0) with driving (£ ^ 0). The bifurcations occur along the contours 
<A 2 ((C 2 - l) 3 ce 6 + 3 (£ 4 + 7£ 2 + 1) a 4 + 3 (£ 2 - l) a 2 + l) = 0. The extra dimension shows that 
an increase in the magnitude of the driving parameter £ means that a stronger coupling between 
the oscillator and two-level system is required to cross the bifurcation. The L z components of the 
fixed points are illustrated in figure [2] and the bifurcations are clearly visible. Also, the £ = line 
correctly reflects the bifurcation shown for the no driving case in figure [I] 



Class 2: for £ 7^ there is one (trivial) fixed point that occurs for all parameter values: 

UJ 2 + K 2 

y° = ^- 

UJ 2 + K 2 



Ll = (21) 
Ll = 



L° z = 

This fixed point is stable for all coupling values. 



3. General qubit dissipation 

The general case of qubit dissipation here means with spontaneous emission present 
(7 > 0, r > 0, k > 0), in which case there are three classes of semi-classical steady states. 
Two of these classes require £ = 0; and the other requires £ ^ 0. 
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Class 1: for £ = there is one fixed point that occur for all parameter values: 



X 




TjU! 


e 




CO 2 + K 2 




y° 




KT] 


K€ 




U 2 + K 2 


AXuj 






1 






= 








= 







(22) 



This fixed point is always stable for /i < v and always unstable for \i > v. 

Class 2: also for £ = there are two fixed points that only occur for a > (3 (note that for 
a — (5 this second class of fixed points is also the first class just described): 



x' = ^ ;r T V2 9 ■ \ -=-^- = -— T V2 



quj r- \uj a — j3 e r- Xuu a — (3 



Lu 2 + k 2 uu 2 + k 2 V Sa 2 4A cj 2 + k 2 V 5a 2 

n 

2/ 



kt\ ^ a — (3 ^2 ftA ja — j3 



lu 2 + ft 2 cj 2 + ft 2 V 5« 2 4Acj cj 2 + ft 2 V 5a 2 

L° = (23) 



L° = ±V2j a - 13 



5a 2 

These fixed points can be stable or unstable depending on the coupling values. 

The L z components of the first two classes of fixed points are plotted in figure [4j The 
bifurcations of these fixed points are shown in the bifurcation diagram of figure [5| 

Class 3: for £ ^ there are up to three real fixed points dependent on a cubic equation: 

id 2 + K 2 UJ 2 + K 2 Z 



UJ 2 + K 2 LO 2 + K, 2 Z 



P L° z (24 ) 



L° = L°, 



where L° z satisfies the cubic equation: 

1 



2^(^-0-^(^-0+^ = (25) 
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Lz((i,v,£=0) with spontaneous emission 



Lz 



1 1 



(i r , 




FIG. 4: L z component of the semi-classical steady states as a function of the parameters \i and 
v when spontaneous emission is present (7 > 0) but there is no driving (£ = 0). There is one 
solution L? z — which is always stable for < v and always unstable for ji > v. The two new 



solutions which appear for ji > v have L z components L® = ±\/2y and can be stable or 
unstable depending on the coupling values. Hence there is a pitchfork bifurcation along the line 
H = v which is often supercritical depending on the coupling parameters. This bifurcation is shown 



Note that for £ 7^ 0, L° z = £ is never a solution to this equation and so the pole in the 
expression for L° above is never encountered. 

If we consider the third class of fixed points at the forbidden point £ = 0, this third class 
of fixed points gives the first (except for L°) and second classes of fixed points; hence it 
generalises the first two classes in a sense. 

The L z component of this third class of fixed points is a function of the three parameters 
/x, z/, and £. The bifurcations of these fixed points are shown in the bifurcation diagram of 



III. QUANTUM STEADY STATES. 

Knowing the coupling parameter values that result in a semi-classical bifurcation of the 
steady state solutions, we wish to investigate whether there is a correspondence with the 
full quantum version. We do this numerically and observe the steady state phase space of 
the oscillator as we change the coupling parameters to move through the semi-classical fixed 




explicitly in the bifurcation diagram of figure [5j 



figure [6) 
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D(//,v,£=0)=0 with spontaneous emission 




1 2 3 4 5 



FIG. 5: Bifurcation diagram of the semi-classical steady states as a function of the parameters 
H and v when spontaneous emission is present (7 > 0) but there is no driving (£ = 0). The 
bifurcations occur along the line \i — v. The extra dimension shows that an increase in the qubit 
dissipation parameter v means that a stronger coupling between the oscillator and two-level system 
is required to cross the bifurcation. The L z components of the fixed points are illustrated in figure 
[4] and the bifurcation is clearly visible. 



D(/i,v,f )=0 with spontaneous emission 




FIG. 6: Bifurcation diagram of the semi-classical steady states as a function of the parameters \i 
and v when spontaneous emission is present (7 > 0) with driving (£ 7^ 0). The bifurcations occur 
along the contours — 2z/£ 4 /x 6 + (/x 2 — 20z//x — 8u 2 ) £ 2 ^ 4 + 8 {ji — z/) 3 fi 2 = 0. The extra dimensions 
show that an increase in either the qubit dissipation parameter v or the magnitude of the driving 
parameter £ means that a stronger coupling between the oscillator and two-level system is required 
to cross the bifurcation. Note that the £ = cross-section correctly reflects the bifurcation shown 
for the no driving case in figure [5j 
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point bifurcation. It is hoped that our semi-classical analysis of the fixed points can be 
numerically justified by observing a signature of the semi-classical bifurcation. 

Here, we use the Quantum optics MATLAB toolbox [10J and pass through two semi- 
classical bifurcations: one by varying the oscillator-qubit coupling (A); and another by 
varying the spontaneous emission (7). By holding all other couplings equal and ignoring 
dephasing, varying these two parameters directly corresponds to varying the parameters 
H and v respectively. Specifically, we will look at a Jahn- Teller model where: u = 0.01, 
A = 0.1, k = 0.001, r = 0, r] = 0, and e = 0. Thus the three parameters on which the 
semi-classical bifurcation depends become: /1 = 3960. 4A 2 , v — 1 + 257 2 , and £ = 0. The 
contour in the bifurcation diagram of figure [5] can thus be redrawn as a function of A and 
7. This is done in figure [7[ 

The MATLAB quantum optics toolbox gives us a steady state density matrix for the 
oscillator-qubit system. From this we can view the steady state phase space of the oscillator 
by plotting the Q-function for the corresponding reduced density operator of the oscillator. 
This is defined [TT] as the matrix elements of the reduced density operator for the oscillator in 
the coherent state basis, Q(a) = tr(p|a)(a|) where \a) is a oscillator coherent state. Three 
series of Q-functions are plotted varying A for two differing fixed values of 7, and varying 7 
for a fixed value of A. These are shown in figures |5J|9j and [TO] respectively. The semi-classical 
bifurcation is clearly evident in each case. 

IV. PHYSICAL IMPLEMENTATIONS OF THE JAHN-TELLER MODEL 
A. A nanomechanical qubit. 

We consider a nonlinear nanomechanical resonator (NR), with resonant frequency uj m) 
coupled to a superconducting microwave resonator, with resonant frequency uo C) see figure [TT| 
The NR is driven parametrically, while the cavity field is driven with two microwave tones 
at frequencies uj\ and uj 2 such that (cj 2 — u\)/2 = uj m and with corresponding amplitudes 
ei, e 2 . The Hamiltonian is 

2 

H = hw c a)a + huo m b ] b + h ^(e*ae^ + e^e"^) (26) 

i=i 

+h X b ]2 b 2 + h(^b 2 e 2lUmt + ^ 2 e~ 2lUmt ) + hg a^a(b + 6 f ) 
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D(A,y)=Ofor fj=0.01,A=0.1,/y=0.001,r=7y=6=0 




7=0.001 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 



FIG. 7: Bifurcation diagram of the semi-classical steady states as a function of the parameters 
A and 7 for the parameters used to numerically calculate the Q-functions: the first series of Q- 
functions in figure [8] is for A increasing along the lower horizontal dotted red line 7 = 0.001; the 
second series of Q-functions in figure [9] is for A increasing along the upper horizontal dotted red 



line 7 = 0.1; and the third series of Q-functions in figure 10 is for 7 increasing along the vertical 
dotted red line A = 0.0225. 




FIG. 8: Q-functions for increasing A along the lower horizontal dotted red line of [7| The sponta- 
neous emission is 7 = 0.0001 giving a qubit dissipation parameter value of v — 1. The semi-classical 
critical value of the oscillator-qubit coupling is A c = 0.01589 at the critical parameter value ji c = 1. 
The steady state phase space of the oscillator is seen to undergo a bifurcation which corresponds 
to the studied semi-classical bifurcation. 
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BBBQD 

FIG. 9: Q- functions for increasing A along the upper horizontal dotted red line of [7} The spon- 
taneous emission is 7 = 0.1 giving a qubit dissipation parameter value of v — 1.25. The semi- 
classical critical value of the oscillator-qubit coupling is A c = 0.017766 at the critical parameter 
value \i c = 1.25. The steady state phase space of the oscillator is seen to undergo a bifurcation 
which corresponds to the studied semi-classical bifurcation. 




FIG. 10: Q-functions for increasing 7 along the vertical dotted red line of Figure [7} The oscillator- 
qubit coupling is A = 0.0225 giving a parameter value of fi = 2.005. The semi-classical critical 
value of spontaneous emission is 7 C = 0.20049 at the critical qubit dissipation parameter value 
v c = 2.005. The steady state phase space of the oscillator is seen to undergo a bifurcation which 
corresponds to the studied semi-classical bifurcation. 
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cavity driving 



co planar cavity 

FIG. 11: Schematic of proposed nanomechanical qubit system. 

where a, a) are the lowering and raising operators for the cavity mode and 6, are the 
lowering and raising operators for the NR. The term proportional to \ represents a quartic 
nonlinearity in the elastic potential energy of the NR and gives rise to a Duffing oscillator [HI 
9j. The term proportional to k represents the parametric driving of the NR and has been 
discussed by [5j . The term proportional to g represents the capacitive coupling between the 
nanomechanical resonator and the cavity field, expanded to linear order in the displacement 
of the nanomechanical resonator. 

We now move to an interaction picture for the NR at frequency cu m and for the microwave 
cavity at frequency (uoi + cj 2 )/2. The Hamiltonian then becomes, 

2 

H = Huj^a + h^2(e*ae iSit + e^V**) (27) 
i=i 

+h X b ]2 b 2 + h(^b 2 + Kb ]2 ) + hg a^a(be- lu;mt + tfe™ mt ) 

where the detuning of the cavity resonance is uj = uj c — (uoi +c^ 2 )/2, and 8\ = (uoi — cj 2 )/2 = 
— 5 2 . Following Woolley et al. [5j we linearise around the steady state amplitudes for 
the cavity field and choose 5i = —uj m) with ei = — e* 2 = ee~ 1 ^ and we find the effective 
Hamiltonian 

H e = hwa^a + hg(a + a t )(6e^ + b ] e 1 ^) + H NR (28) 

where 

H NR = h X b ]2 b 2 + hK{b 2 + b ]2 ) (29) 
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with k real, and g = goe/uj m) and where \± is the decay rate for the cavity field. Note that 
the coupling constant, g, can be made large by increasing the driving field e. 

The dynamics arising from H NR was investigated in [7j. The classical phase space dy- 
namics has two elliptic fixed points either side of a hyperbolic unstable fixed point. This is 
equivalent to an effective double well potential. The quantum ground state is then seen to 
be very well approximated by a symmetric superposition of two oscillator coherent states 
| ± a) centered on the fixed points with 

a = -<y| (30) 

The energy separation between the ground state and first excited state was shown in jTj to 
be 

AE = h X e~ 2lal2 = HA S (31) 

We now assume that the NR is always very close to its ground state and we truncate the 
Hilbert space to the ground state and the first excited state. We then define a qubit basis 
by 

|0> = \a) (32) 
|1> = | - a) 

It is then clear that 6|0) = a|0), 6|1) = — a\l). If we then define a z — — |0)(0|, a x = 

|1)(0| + |0)(1|, we can then write 

H NR h—a x (33) 
and the effective Hamiltonian, with the phase choice ifj = 7r/2, takes the form 

H e = hjjo^a + h^-a x + h\(a + a))(j z (34) 

where 

K + ^2) 

uj = uj c (35) 

A = 2g\a\ (36) 

A = 2 X e" 2|a|2 (37) 

This is the Jahn- Teller model with a critical coupling strength X cr given by[6j 

X cr — - — (38) 
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If we add a resonant driving term to the NR, we get an additional term proportional to a z . 

In [8], the elastic nonlinearity of a Pt nanowire, implemented by Kozinsky et al.[9j, was 
estimated to have \ = 10 -4 s -1 . If we keep the parametric driving weak, we can ensure that 
|a| = 5 (see [5j), this is a deep quantum regime in which the fixed points are separated from 
the unstable fixed point at the origin by a few units of ground state uncertainty. If we accept 
a detuning of about 10MHz, uj ~ 10 8 , then X cr ~ 0.7s -1 . The coupling constant is largely 
determined by g which depends on the intra-cavity mean field amplitude. This quantity 
can thus be controlled quite well. Typical values [5] for g are of the order of 10 — 1000s -1 , 
with the smaller number for bad cavities, so it would be relatively easy to exceed the critical 
coupling strength. 



B. Circuit QED 



Devoret et al. [T] have proposed a scheme to get ultra strong coupling between a Cooper 
pair box qubit and the microwave field of a coplanar superconducting resonator. The central 
conductor of the coplanar cavity is divided into two segments separated by a Cooper pair 



box, see figure [12] The quantum theory of such a system begins by first writing down 
the classical circuit dynamics, constructing a Lagrangian and an assocaited Hamiltonian. 
Quantistation then proceeds via the usual canonical method. This results in an effective 
quantum theory in which collective variables of direct interest to the experimentalist couple 
only weakly to the microscopic degrees of freedom, which remain as a source of dissipation 
and decoherence. 




FIG. 12: Lumped element schematic circuit approximation. 
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We label the relevant electrical variables as shown in figure [T3J. We consider an input 
and output currents ii(t) and i 2 (t) at voltages v\{t) and v 2 (t) respectively. The "mirror" 
capacitors of capacitance Cmi and Cm2 hold charges Qmi(*) an d Qm2(*) respectively and 
enable the device to be inserted into a transmission line. The lumped inductances ^ of 
the "cavity" hold magnetic fluxes $i(i) and $2^)- The lumped capacitances ^ of the 
cavity hold charges Qi(t) and Q2(t). The pure Josephson element, with a critical current 
I j and tunneling energy Ej, has a wavefunction phase difference across it of 6{t) and sees 
N(t) Cooper pairs tunnel across it. The Josephson junction capacitance Cj holds a charge 
Qjc{t). The currents flowing through the inductive and capacitive elements of the cavity are 
iL(t) and ici(t) and ic2(t) respectively. The currents flowing through the pure Josephson 
element and the Josephson junction capacitance are ij(t) are ijc(t) respectively. Finally, 
the voltage at the input and output ends of the resonator are vm(t) and VR2(t) respectively; 
and the voltage at the input and output ends of the Josephson junction are vji(t) and Vj2(t) 



respectively. Figure 13 summarises all of this. 
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FIG. 13: Lumped element schematic circuit approximation, here labelled with the dynamic elec- 
trical variables. 

If we introduce creation and annihilation operators of the cavity field (the effective flux 
$ acts as a position, and the effective charge Q as a momentum): 
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where: 

1 1 1 

+ 



Ce C + Cm Cj 
1 



(40) 



\JLCe 

and if we consider only the bottom two energy levels for the Josephson junction using Pauli 
matrices for the resulting qubit, then we can rewrite the Hamiltonian (ignoring constant 
energy offsets) in the Jahn- Teller form of Eq.([3]) where, following Devoret et al.pQ, 

2A If E c \ 1/4 /Z vac _ 1/2 



(42) 



A ^ \2EjJ V Z c 

with 

Z va c = — ~ 377f2 

Typically Zc = 50 Q. Taking Ec/Ej = 200, Devoret et al.[lj arrive at a value of A w 10 A 
which is certainly well outside of the domain of validity for the rotating wave approximation. 
We thus believe this configuration offers a good chance of designing a system with a coupling 
strength above the Jahn- Teller dissispative bifurcation in circuit QED. 



V. CONCLUSION. 



In this paper we have presented a detailed analysis of the effect of dissipation on the 
dynamical bifurcation that occurs when there is strong coupling between a single two level 
system and an oscillator; a Jahn- Teller model. We have based our description of dissipation 
on the physically appropriate mechanisms two physical realisations of the system based on 
nanomechanicas and circuit Quantum electrodynamics. The key feature of the bifurcation 
in the dissipative Jahn- Teller model is the change in the oscillator fixed point from one 
centered on a point of zero radius in phase-space to one with support on a non zero value of 
the radius. This is a distinct kind of bifurcation from that discussed recently in the damped 
nanomechanical Duffing oscillator [9] which only involved a single degree of freedom. In the 
case considered here, the bifurcation results in steady state correlations between the state 
of the oscillaator and the two-level system. 

As the average excitation energy an oscillator is proportional to the radius in phase-space, 
this bifurcation would be reflected in a change in the steady state mean excitation energy 
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from zero to a finite non zero value. This would have implications for any attempt to cool 
the system though tuning to the red sideband transition, that is to say, tuning the cavity 
field driving by the mechanical frequency below the cavity frequency. If the parameters were 
such that the system was already beyond the Jahn- Teller bifurcation, the mechanical system 
could not be cooled to a zero phonon state, but would rather relax to the bistable state with 
a non zero mean phonon number. Fluctuations would then drive switching events between 
the two stable steady states. 

The non-dissipative model, for coupling stronger than the critical coupling, has a ground 
state with significant entanglement between the two-level system and the oscillator. We do 
not know if any entanglement remains in the steady state of the dissipative model beyond 
the bifurcation point. This is a difficult question to answer as the steady state has a non 
Gaussian Q-function (or Wigner function) and thus it is not clear what would be a good 
measure of entanglement. In a future work we will use positive P-function methods to 
attempt to answer this question. 

In many implementations of quantum information processing there is often an unwanted 
strong coupling between an oscillator degree of freedom and a strongly damped two-level 
system [131 CH]- The model of this paper may be relevant to the on going study of such 
systems in those cases where a perturbative treatment of the coupling is not possible. 

We have given two examples of quantum electromechanical systems that could exhibit the 
steady state bifurcation of the dissipative Jahn- Teller model. Observation of this effect would 
be a clear demonstration of the ultra strong coupling regime that can be achieved in these 
systems, as opposed to what typically happens in atomic systems where the rotating wave 
approximation eliminates the bifurcation. Such system open a path to study the quantum 
signature of non linear bifurcations on the steady states of strongly coupled systems. 
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